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We present a measurement of the branching fraction for the decay B~ — > D°K*~ using a sample 
of approximately 86 million BB pairs collected by the BABAR detector from e + e~ collisions near the 
T(4S) resonance. The D° is detected through its decays to K~ir + , K~-k + iv and K~tt + iv~it + , and 
the K*~ through its decay to K% tt~ . We measure the branching fraction to be B(B~ — > D°K*~) 
= (6.3±0.7(stat.)±0.5(syst.))xl0^ 4 . 
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PACS numbers: 13.25.Hw 14.40.Nd 



A comprehensive test of CP violation within the 
Standard Model requires precision measurements of the 
three sides and three angles of the Unitarity Triangle, 
which are combinations of various Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements Q. The measurement 
of the angle 7 of the Unitarity Triangle is challenging and 
requires larger samples of B mesons than are currently 
available. A precise determination of 7 at the B facto- 
ries is likely to use many different decay modes. Decays 
of the form B -> D^K^ can provide a theoretically 
clean determination of 7 ■ For some of the proposed 
methods, there are distinct advantages to using the K* 
modes 3]. In this paper, we measure the branching frac- 
tion for one of these decays, B~ — > D°K*~ which was 
first observed by the CLEO experiment [5j- If the D° is 
reconstructed in its decay to CP eigenstates, the b — > cus 
and b — > ucs quark transitions interfere, giving access to 
the phase 7 through the measurement of direct CP vi- 
olation asymmetries. However, the branching fractions 
for D° decays to CP eigenstates are only of the order of 
1%, too small for the size of the available data sample. 
Therefore, for this analysis, we use decay modes of the 
D° and K*~ that have clear experimental signatures and 
sufficiently high branching fractions. This measurement 
provides an important first step towards establishing the 
feasibility of using the decay B~ — > D°K*~ for a future 
determination of 7. 

We present here a measurement of the branching frac- 
tion for the decay B~ — > D°K*~ using data collected 
with the BABAR detector at the PEP-II e + e~ storage 
ring. The data correspond to an integrated luminos- 
ity of 81.5 fb _1 taken at center-of-mass energies close 
to the T(AS) resonance, giving a sample of approxi- 
mately 86 million BB pairs. We reconstruct D° candi- 
dates through the decays D° -> K~tt + , D° -> X^tt+tt 
and D° — > K~Tr + ir~ir + . K*~ candidates are identified 
through the decay K*~ — > K® tt~ , with the decaying 
to a pair of charged pions. 

A detailed description of the BABAR detector can be 
found elsewhere 6J. Only detector components relevant 
to this analysis are described here. Charged-particle tra- 
jectories are measured by a five-layer double-sided sil- 
icon vertex tracker (SVT) and a 40-layer drift cham- 
ber (DCH), operating in the field of a 1.5-T solenoid. 
Charged-particle identification is achieved by combining 
measurements of ionization energy loss (dE/dx) in the 
DCH and SVT with information from a detector of inter- 
nally reflected Cherenkov light (DIRC). Photons are de- 
tected in a CsI(Tl) electromagnetic calorimeter (EMC). 
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We set the event-selection criteria to minimize the sta- 
tistical error on the branching fraction, using simulations 
of the signal and background. In general, charged tracks 
are required to have at least 12 DCH hits and a mini- 
mum transverse momentum of 0.1 GeV, and to originate 
from the interaction point, within 10 cm along the beam 
direction and 1.5 cm in the transverse plane. We use less 
restrictive selection criteria for tracks used to reconstruct 
K° — ► 7r + 7r~ candidates, to allow for displaced K® decay 
vertices. Photon candidates are identified in the EMC 
as deposits of energy isolated from charged tracks. They 
are required to have a minimum energy of 30 MeV and a 
shower shape consistent with that of a photon. 

We use pairs of photons to reconstruct ir° candidates, 
which are required to have an invariant mass between 125 
and 144 MeV. We reconstruct K® candidates from pairs 
of oppositely charged tracks fitted to a common vertex. 
They are required to have an invariant mass within 8 MeV 
of the K° mass Q. 

To reconstruct K*~ candidates, we combine candi- 
dates with charged tracks. We require the K*~ candidate 
to have an invariant mass within 75 MeV of 892 MeV. In 
addition, the K® vertex is required to be displaced by at 
least 3 mm from the K*~ vertex. 

We reconstruct D° candidates from the appropriate 
combination of tracks and tt° candidates. The K~ tracks 
must satisfy kaon identification criteria resulting in an 
efficiency of 80%-95% depending on the momentum. The 
probability of a pion to be misidentified as a kaon is less 
than 5%. We require the momenta of the K _ candidates 
to be greater than 250 MeV and their polar angle (relative 
to magnetic-field axis) to be in the interval 0.25 < 8 < 
2.55 rad to restrict them to a fiducial region where the 
kaon identification performance can be determined with 
small uncertainty. The tracks from the D° are fitted to a 
common vertex and we accept candidates if they have an 
invariant mass within 18 (14) MeV of the D° mass for 
the K~it + (K~it + n~it + ) decay. For the K~ir + ir° decay, 
we use an asymmetric mass requirement —29 < (m — 
1865 MeV) < +24 MeV, reflecting the distribution of the 
energy of the photons from the tt° decay. It is known that 
the decay D — > A~7r + 7r° occurs predominantly through 
an intermediate state (K*~ (892) or p+(770)). Hence, 
to reduce the combinatorial background in the A~7r + 7r° 
decay, we select events in the enhanced regions of the 
Dalitz plot, using amplitudes and phases determined by 
the CLEO experiment Q- 

In reconstructing the decay chain, the measured mo- 
mentum vector of each intermediate particle is deter- 
mined by refitting the momenta of its decay products, 
constraining the mass to the nominal mass of the parti- 
cle and requiring the decay products to originate from a 
common point. For the K*~ resonance only a geomet- 
rical constraint is used in this kinematic fit. Finally, to 
reconstruct B~ decays, D° candidates are combined with 
A'*~ candidates. 
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The dominant background is from e + e~ — > produc- 
tion. We suppress this background using requirements 
on the event topology and kinematics, and through the 
use of a Fisher discriminant. The ratio of the second and 
zeroth Fox- Wolfram moments 0, which is a measure of 
the event sphericity and is close to zero for approximately 
spherical events, is required to be less than 0.5. The ab- 
solute value of the cosine of the angle between the thrust 
axis of the B candidate and the thrust axis of the rest 
of the event, |cos#t|, is peaked at one for continuum 
events and is approximately flat for B decays. We re- 
quire | cos#t| < 0.8 for K~n + and A"~7r + 7r° decays and 
Icos^tI < 0.75 for K~ir + TT~ir + decays. The Fisher dis- 
criminant is built from the momentum of all particles 
in the event (excluding those used to form the B can- 
didate) and the angle between this momentum and the 
thrust axis of the reconstructed B~, both in the center- 
of-mass frame The K*~ helicity angle, 9h, defined 
as the angle between the tt~ from the K*~ decay and 
the B~ flight direction in the rest frame of the K*~, 
follows a cos 2 9h distribution for signal events and is ap- 
proximately flat for continuum events. To further reject 
continuum background in the K~ir + n~Tr + channel, we 
require |cos##| > 0.4. 

The selection criteria just described reject all but ap- 
proximately 0.001% of the background, while retaining 
between 4% and 13% of the signal, depending on the D° 
mode. In the case of events with more than one B~ can- 
didate (5-17%, depending on the D° mode), we choose 
the best candidate on the basis of the \ 2 formed from 
the differences of the measured and true B~ , D°, and 
K® masses, scaled by the mass resolutions. Studies of 
simulated signal events have determined that the algo- 
rithm does not introduce a bias and chooses the correct 
B~ candidate in approximately 80% of the events with 
multiple candidates. 

Finally, we identify _B-meson decays kinematically 
using two nearly independent variables: the energy- 
substituted B mass Toes 



VGs^ + po-Pb) 2 /^ 2 
where the subscripts and B refer to the e + e system 
and the B candidate respectively, s is the square of the 
center-of-mass energy, and energies {E) and momentum 
vectors (p) are computed in the laboratory frame; and 
AE = E* B - \fs/2, where E* B is the B candidate en- 
ergy in the center-of-mass frame. We select B~ candi- 
dates with |A_E| < 25MeV, which corresponds to ap- 
proximately ±2.2cr (where the resolution a is found to 
be independent of the D° decay mode). In addition, the 
signal events are expected to have values of toes close to 
the B~ mass. 

We determine the signal yield of B~ — > D°K*~ events 
by performing an unbinned maximum likelihood fit to the 
toes distribution of the selected candidates for the signal 
region in AE. The signal distribution is parameterized 
as a Gaussian function and the combinatorial background 
as a threshold function All parameters except the 

endpoint of the threshold function are unconstrained in 
the fit. 



The signal yield determined from the fit potentially 
includes backgrounds from other BB decays that also 
peak in toes- To investigate this, we have studied a sim- 
ulated sample of generic BB decays and also high statis- 
tics simulated samples of other B — > d(*)AT(*) decays. 
The simulation indicates no enhancement in the signal 
region from this background. Therefore, we assume that 
the peaking background is negligible and the uncertainty 
in its determination from the studies of various simulated 
event samples is included as a systematic error. We have 
also verified that use of the B~ mass and error in the \ 2 
calculation for the choice of the best B~ candidate does 
not affect the smooth shape of the background in toes- 

Figure ^ shows the Toes distribution for the three dif- 
ferent D° decay modes with the fit function superim- 
posed. A clear signal is seen in all cases. The signal 
yield and the size of the combinatorial background in 
the signal region, i.e., for events with toes > 5.270 GeV, 
are detailed in Table [I] We observe a total of 161 ± 17 
B~ — > D°K*~ events. We have studied the cos 9h distri- 
bution for the selected candidates and determined that 
the data are consistent with pure B~ — > D°K*~ decay. 
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FIG. 1: The m E s distributions of B~ -> D°K*~ candi- 
dates: (a) D° -> R-TT+, (b) D° -> R-tt+tv , and (c) 
D° — » A' - 7r + 7r~7r + . The solid lines show the fit used to ex- 
tract the signal yields, with the distribution parameterized 
as a Gaussian plus a threshold function as described in the 
text. The dashed line indicates the combinatorial background 
component. 
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TABLE I: Signal yield, number of background events, and 
efficiency for the three D° decay modes in the signal region 
(m,ES > 5.270 GeV). Yields are extracted from the fits to 
the rriES distribution from data (errors are statistical only). 
Efficiencies are computed from simulated events. 





K'TY + 


K-TV + TT { > 


K 7T + 7T 


-7T+ 


Signal Yield 


56.2 ± 9.4 


51.7 ± 11.0 


52.6 ± 


8.7 


Background 


f9.5 ± 4.3 


37.7 ± 6.2 


16.4 ± 


3.6 


Efficiency (%) 


f2.8 


3.5 


4.0 





We determine the selection efficiency for each sample 
of B~ — > D a K*~ events from samples of simulated signal 
events. We apply small corrections determined from data 
to the efficiency calculation to account for the overesti- 
mation of the tracking and particle-identification perfor- 
mance, and of the tt° and K® reconstruction efficiencies 
in the Monte Carlo simulation. The product of these 
efficiency corrections is about 0.9. 

To quantify the ability of the simulation to model 
the variables used in the event selection, we use a sam- 
ple of B~ — > D°tt~ events from data and Monte Carlo 
simulation. This sample is kincmatically similar to the 
B~ — > D°K*~ decay. We select B~ — > D°tt~ events in 
the same way as the B" — > D°K*~ sample, with the 
additional requirement that the 7r~ fails loose kaon iden- 
tification criteria, to remove B~ — > D°K~ events. Ap- 
proximately 3000 B~ — > D Q ir~ candidates in each D° 
decay mode are selected from the data. The purity of 
the sample is 94% for the K~ir + ir° decay and 98% for 
the K~tt + and K~Tr + ir~Tr + decays. We use this sample 
to determine correction factors for the efficiencies for the 
B~ — > D°K*~ selection. The obtained correction factors 
vary from about 0.95 for the K~ir + and K~Tr + ir~ir + de- 
cays to 0.85 for the K ~7r + 7r° decay. We include the sta- 
tistical precision of these corrections in the systematic 
error of the branching fraction. The selection efficiency 
after all corrections are reported in Tabled 

We determine the branching fraction separately for 
each of the D° decay modes from: 



for a signal yield of N events, an efficiency e and a sample 
containing N BB - pairs of B mesons. Spo, Bjc*- > ^k° s an d 
B^o, the branching fractions for the D°, K*~ , K® and 
7T° respectively to the relevant final states, are obtained 
from Ref. ,,.0 in equation is only relevant for the 

K~ir + ir° mode). We assume that the T(AS) decays to 
pairs of B + B~ and B°B° mesons with equal probability 
and we do not include any additional uncertainty due to 
this assumption. 

We have identified several sources of systematic un- 
certainty as significant. The number of BB pairs in the 
data sample is known with an uncertainty of 1.1%. The 
uncertainties in the D° branching ratios are taken from 



Ref. |3 ■ We determine the systematic errors arising from 
uncertainties in track, K® and ir° reconstruction and in 
kaon identification from studies of high statistics data 
control samples. The uncertainty in the track recon- 
struction efficiency is determined to be 0.8% per track 
originating from the interaction region. There is an ad- 
ditional uncertainty of 3% arising from the knowledge of 
the K® reconstruction efficiency. The charged kaon iden- 

TABLE II: Systematic uncertainty estimates for each of the 
three D° decay samples. 



Source Uncertainty (%) 









K 7r + 7r 7r + 


Number of BB events 


1.1 


1.1 


1.1 


Simulation statistics 


•5.6 


6.5 


6.0 


D° branching ratio 


2.4 


6.2 


4.2 


Tracking efficiency 


2.4 


2.4 


4.0 


Ks efficiency 


3.0 


3.0 


3.0 


Particle identification 


2.0 


2.0 


2.0 


7T° efficiency 




5.0 




Peaking background 


2.3 


1.4 


3.1 


K*~ lineshape 


3.0 


3.0 


3.0 


Data/simulation differences 


1.4 


2.4 


2.1 


Total 


8.6 


11.9 


10.3 



tification leads to a systematic uncertainty of 2%, and 
the 7T° reconstruction to a systematic uncertainty of 5%. 
The systematic error from the knowledge of the peaking 
background is taken from the studies of various simulated 
data samples described above. An additional uncertainty 
from the knowledge of the K*~ lineshape has been de- 
termined to be 3%. Finally, we include the errors on the 
correction factors determined from the B~ — > D°ir~ sam- 
ple. We have studied the uncertainty in the parameteri- 
zation of the background and of the signal by repeating 
the 771es fits with different combinations of parameters of 
the functional form fixed to values obtained either from 
simulation or from studies of sideband regions in AE. 
We conclude that the systematic uncertainty from this 
source is negligible. A summary of the systematic errors 
estimate is shown in Table ITT1 

The resulting B branching fractions corresponding to 
three different D° decay modes are listed in Table Mil 
We determine the weighted average of the three mea- 

TABLE III: Measured branching fraction B(B~ -> D°K*~). 
The first errors are statistical and the second systematic. 



Decay Mode g(10~ 4 ) 
D u ->• A'~7r+ 5.8±1.0±0.5 
D° -» K'-k+tt 5.8±1.2±0.7 
D° -> A'~7r+7r~7r + 8.7±1.5±0.9 
Weighted Average 6.3±0.7±0.5 



surements, B(B~ -> D°K*~) = (6.3 ± 0.7 ± 0.5) x 10~ 4 , 
taking into account the correlations between the sys- 
tematic uncertainties. The result of this analysis is in 
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good agreement with a previous measurement by CLEO, 
B(B- -> D°K*-) = (6.1 ± 1.6 ± 1.7) x 10" 4 @. 

In summary, we have studied the decay B~ — » D°K*~ , 
where the D° was detected through its decays to K~ir + , 
K~ir + ir° and K~ir + Tr~ir + and the K*~ through its de- 
cay to K ° tt~ . We have measured the branching fraction 
B{B- — > D°K*~) = (6.3 ± 0.7 ± 0.5) x lO" 4 . This is in 
good agreement with the previous measurement of this 
branching fraction, and significantly improves on its pre- 
cision. In the future, with larger data samples, this decay 
will be studied with the D° reconstructed in CP eigen- 
states. Eventually it is hoped that decays of the form 
B — ► D^*'K^ can provide important constraints on the 
angle 7 of the Unitarity Triangle. 
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